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ABSTRACT 

Hadronic emission from supernova remnant (SNR)-molecular cloud (MC) as¬ 
sociation systems has been widely regarded as a probe of the shock-accelerated 
cosmic-ray protons. Here, we report on the detection of a 7 -ray emission source 
with a significance of 24a in 0.2-300 GeV, projected to lie to the northwest of 
the thermal composite SNR Kesteven 41, using 5.6 years of Fermi -Large Area 
Telescope (LAT) observation data. No significant long-term variability in the en¬ 
ergy range 0.2-300 GeV is detected around this source. The 3a error circle, 0?09 
in radius, covers the 1720MHz OH maser and is essentially consistent with the 
location of the Wsr ~ —50 km s ' 1 MC with which the SNR interacts. The source 
emission has an exponential cutoff power-law spectrum with a photon index of 
1.9±0.1 and a cutoff energy of 4.0±0.9 GeV, and the corresponding 0.2-300 GeV 
luminosity is ~ 1.3 x 10 36 ergs _1 at a distance of 12kpc. There is no radio pul¬ 
sar in the 3a circle responsible for the high 7 -ray luminosity. While the inverse 
Compton scattering scenario would lead to a difficulty in the electron energy bud¬ 
get, the source emission can naturally be explained by the hadronic interaction 
between the relativistic protons accelerated by the shock of SNR Kesteven 41 
and the adjacent northwestern MC. In this paper, we present a list of Galac¬ 
tic thermal composite SNRs detected at GeV 7 -ray energies by Fermz-LAT is 
presented. 
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Subject headings: acceleration of particles — gamma rays: ISM — ISM: super¬ 
nova remnants — ISM: individual objects (Kes 41) 


INTRODUCTION 


The origin of the cosmic rays (CRs) remains one of the most debated issues in high 
energy astrophysics more than 100 years after they were detected for the first time. Super¬ 
nova remnants (SNRs), whose strong shocks contain huge amounts of energy, are consid¬ 
ered to be t he most pro bable candidates among the Galactic CR acceleration sources (e.g., 
Ginzburg fc Svrovatsk ii 1969|). A long standing argument concerning the putative SNR-CR 
link is that non-thermal radio emission from SNRs provides clear evidence for electron ac¬ 
celeration, whereas the CR spectrum observed on Earth is 99% protons and other nuclei. 
Therefore, the y-ray emission from SNRs that is dominated by the decay of n° mesons 
produced via proton-proton collisions (i.e., the hadro nic in tera ction) plays a key role in pro¬ 
viding evidence for proton acceleration flAckermann et al.ll2013l) . However, it is often difficult 
to distinguish between the hadronic y-ray emission and the electrons’ inverse Compton or 
non-thermal bremstrahlung emission (i.e., the leptonic emission). There are generally two 
scenarios that describe how hadronic y-rays are produced in SNRs. In one scenario, the 7r°- 
decay emission is suggested to arise from shock-crushed dense clouds where th e accelerated 
protons frozen in the clouds efficiently collide with ta rget cloud gas (e.g.. iBlandford fc Cowie 
19821: lUchivama et al.ll2010l: iTang fe Chevalierl 120141) . In the other scenario, the hadronic y- 
rays are ascribed to interactions between th e relativistic protons esca ping f r om the SNR 


shoc k and adjacent m olec ular clouds (MCs; lAharonian fe Atovanl e.g., Il996l: iGabici et al. 


e.g. 


20091 : iLi fc Chenl e.g., l2010l: lOhira et al 


e-g., 


20111). In both scenarios, the SNRs inter¬ 


acting with MCs are crucial probes in the search for the signatures of proton acceleration. 
The hadronic y-ray emission from SNR-MC systems is usually bright around GeV, and a se¬ 
ries of GeV-bright SNRs interacting with MCs have recently been discovered with the Large 
Area Telesco pe (LAT) on boa rd the Fermi Gamma-ray Space T elescope. These SN Rs in- 
clude W51C ( Abdo et al. 20091). W44 (Abdo et al. 2010c ), IC 443 (Abdo et al. 2010d ). W28 


f Abdo et al. 2010a ). W41 f Castro et al. 20131 ). RCW 103 ( Xing et al 


2 0141) . etc. Additional 


GeV observations continue to enlarge the sample of hadronic interaction between SNRs and 
MCs, and here we present a GeV study of another SNR, namely, Kesteven 41(G337.8—0.1). 

As a southern-sky SNR, Kes 41 is shown to be centrally brigh tened in X-rays within a dis¬ 
torted radio shell by a XMM-Newton observation flCombi et al.ll2008l). and th e refore is clas¬ 
sified as a thermal composite (or mixed-morphology) SNR fl Jones et al.lll998l : [Rho & Petre 
19981 ). The X-ray emitting plasma of the SNR has been newly revealed to be rich in sulfur 
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and argon; thus, Kes 41 jo ins the subcl ass of “enhanced-abundance” or “ejecta-dominated” 
thermal composites (j Zhang et al . 2 0151 ). Kes 41 has also been found to be interacting with 
an adjacent MC, as indicated by the 1720 MHz hydroxyl radical (OH ) maser emission de¬ 
tected in the northern radio shell flKoraleskv et al.l [19981 ; ICaswel l 20041). Recently, we found 
that Kes 41 is associated with a giant MC at a systemic local standard of rest (LSR) ve¬ 
locity of —50kms _1 and is confined in a cavity deli neated by a northern molecular shell, a 
western concave MC, and a southeastern H I cloud flZhang et al.ll2015l ). The forward shock 
is suggested to have left the adiabatic stage since the SNR shock encountered the cavity 
wall, while the inner thermal X-rays are ascribed to heating by the reflection shock from 
the cavity wall. The birth of Kes 41 inside the molecular cavity provides a mass estimate of 
> 18M 0 for the stellar progenitor. It is logical and meaningful to search for the hadronic 
emission due to the interaction of the SNR with the dense environmental gas. 


In this paper, we report the results from a spatial and spectral analysis of the Fermi- 
LAT observation data of the Kes 41 region. We describe the Fermi observation data in 
Section [2] and present the data analysis and results in Section [3j The possible physical 
relation of the detected y-ray emission with the SNR is discussed in Section [JJ 


2. OBSERVATIONS AND DATA REDUCTION 


The LAT on board Fermi, launched on 2008 June 11, is a y-ray imaging instrument that 
covers a very wide range of energy from 20 MeV and up to 300 GeV. It reconstructs the direc¬ 
tion of incident y-rays by tracking the electrons and positrons resulting from pair conversion 
of the y-rays in the solid state silicon trackers, and measures the ener gy o f th e su bs equen t 
electromagnetic showers that develop in the cesium iodide calorimeters (lAtwood et al.l 20 091). 
The point-spread function (PSF) varies largely with photon energy and im proves a t high en¬ 
ergies (the 68% containment radius at >2 GeV is smaller than 0?5 [Atwood et al. 20091 ). 


We use the reconstructed Pass 7 reprocessed version of 5.6 years of accumulated Fermi- 
LAT data@ that has been selected from 2008 August 04 15:43:37 (UTC) to 2014 April 
01 02:29:28 (UTC). We analyze the data with the standard software, ScienceTools ver¬ 
sion v9r32pfo released on 2013 October 24, with the instrument response functions (IRFs) 
P7REP_SOURCE_V15. Standard selection criteria are applied to the data selection process 
as described below. The Source (evclass=2) events are selected and the maximum zenith 


1 http://fermi.gsfc.nasa.gov/ssc/data 
2 See http://fermi.gsfc.nasa.gov/ssc 
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angle cut is 100° to reduce the residual 7 -rays from CR interactions in the upper atmosphere. 
We used the standard criteria for selecting time intervals for analysis: (DATA_QUAL==1) 
&& (LAT_CONFIG==l) && ABS(ROCK_ANGLE)<52. The analysis is restricted to the 
energy range above 200 MeV due to uncertainties in the effective area and broad PSF at low 
energies and below 300 GeV due to limited statistics. 


ANALYSIS AND RESULTS 


In our analysis, we select the LAT events inside a 14° x 14° region of interest (ROI, 
in equatorial coordinate system) centered at the position of the Kes 41 (R.A. (J2000) = 
16 h 39 m 00 s and decl. (J2000) = —46°58 , 59") with a bin size of 0°.04 x 0°.04. We per¬ 
form our analysis following the stand ard b inned likelihood analysis procedure. The second 
Fermi -LAT Catalog (2FGL) sources flNolan et al. 2Q12) within radius 15° around Kes 41 
are included in the source model, which was generated by the user-contributed software 
make2FGLxml. pj|f]. The Galactic and extragalactic diffuse background components (as spec¬ 
ified in the hies glLiem-v05.fits and isosource05.txt, respectively ) are used. In the likelihood 
fittings, the spectral parameters of the sources located beyond 10° of the ROI center are 
fixed to the values reported in 2FGL, and the spectral parameters of all the sources located 
within 10° of the center of ROI, together with the normalizations of the two diffuse back¬ 
grounds, are allowed to vary. The fittings are performed with the optimizer NEWMINUIT 
until convergence is achieved. 


3.1. Source Detection 

First, a binned likelihood analysis is applied in the energy range 2-300 GeV. In the source 
model, the source 2FGL J1638.0—4703c, which is very close to Kes 41, has been removed due 


to the uncertainty of its spatial and spectral information caused by the imperfectly modeled 

diffuse emission, and thus needs to be treated with great cared ( 

Nolan et al 

2012 

). A newly 

discovered 7 -ray source (HESS J1641—463, 

Lemoine-Goumard et al. 

2014) 

has been added 


assuming a power-law spectrum. Then, the test statistic (TS, defined as 2(log£ — log£ 0 )> 
here jC 0 is the likelihood of null hypothesis and C is the likelihood with the source included) 
map for a 1° x 1° region centered at Kes 41 is made after subtracting this baseline model (see 


3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 

4 http: / /fermi. gsfc. nasa.gov/ssc/data/access / lat / 2yr jcatalog / 
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Fig. [[]). As can be seen in Figured! there is excess 7 -ray emission in the region of Kes 41. 
The position of the peak of the TS value is on the northwest of the SNR, but does not agree 
with the position of 2FGL J1638.0—4703c. It is generally consistent with the lo cati on of the 
d ense MC at Wsr ~ — 50kins -1 which is found to be associated with the SNR ( Zhang et ah 


20151 ). We also perform an analysis in the low energy range 0.2-2 GeV, and some residual 


7 -ray emission is detected at the same position. Therefore, we add a point source with a 
power-law spectrum at the position where the TS value is highest in our source model to 
approximate the excess emission. After that, we conduct a binned likelihood analysis in the 
broad energy range 0.2-300 GeV and, utilizing gtfindsrc (a tool in the LAT software package 
ScienceTools ), we find the best-fit position of the excess 7 -ray emission at (R.A. (J2000) = 
16 h 38 m 36100, deck (J2000) = —46° 55' 06".96) with la nominal uncertainty of 0?03 and 
3a nominal uncertainty of 0?09. 


By comparison, we detect this source as a point-like source with a power-law spectrum in 
0.2-300 GeV with a significance of 24a at the best-fit position and increases the significance 
by la over the position of 2FGL J1638.0—4703c. The data we use are collected from 5.6 years 
of Fermi -LAT observations while the tentative source 2FGL J1638.0—4703c was suggested 
based on the first two years of observations. Both the statistical result and increased exposure 
time hence suggest that the 7 -ray emission excess at the best-fit position adjacent to Kes 41 
is more significant than 2FGL J1638.0—4703c. Thus, we replace 2FGL J1638.0—4703c with 
this new source at the best-fit position (hereafter source A) in the following analysis. 


In an attempt to explore the origin of tlm 7 -ray emission of source A, we searched in 
the SIMBAD Astronomical Database (IWenger et al.l 200C8 ) within a 3a error circle of the 
source (see Fig. [[]). In addition to SNR Kes 41, only nine dark clouds, a young stellar object 
candidate, and an infrared source are known to exist in the region. Therefore, the origin of 
this 7 -ray emission is most likely related to the SNR. 


3.2. Timing Analysis 


We next search for long-term variability in the one month binned light curve of sou rce A 
in the energy range 0.2-300 GeV, which is obtained from likelihood analysis (INolan et al. 
2 01 21 ) in each time bin. As can be seen in the light curve (Figure [ 2 ]), all of the flux points 
remain within 3a uncertainties of the average flux. Fitting the flux points with TS value > 4 
to a constant flux model (shown as a red line in Figure [2]) yields a y 2 ~ 32.3 with 48 degrees 
of freedom (dof). Moreover, we calculate the Variability Judex, TS var , of source A (with all 
69 time bins) in the 0.2-300 GeV energy range according to the method introduced in §3.6 
of Nolan et ah ( 20121 ). If the flux is constant, then TS var is distributed as y 2 with 68 dof, 
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and variability would be considered probable when TS var could exceed the threshold of 98.0 
corresponding to 99% confidence. The computed TS var of source A is 65.2, corresponding to a 
confidence level < 50% for a variable source. These results suggest that there is no signihant 
long-term variability observed in the region of Source A in the 0.2-300 GeV energy range. On 
the other hand, we construct 1000 s binned light curves of source A in the same energy range 
which are obtained through Fermi-LAT aperture photometry analysed using LAT photons 
within different aperture radius from 0?2 to 0?5 . We analyze these light curves for periodic 
signals, but no significant periodicity is detected. However, this method is statistically 
limited and the periodicity is hard to detect due to the massive diffuse background photons 
in a low galactic latitude. 


There is a close positional correspondence between source A, suggested here as a steady 
source, and 3FGL J163 8.6—4654, w hich is indicated as a variable source in the third Fermi- 
LAT Catalog (3FGL) f Acero et al. 2015 1. 3FGL J1638.6—4654 is detected in 0.1-300 GeV 
with a significance of 13cr and a 0.1-300 GeV energy flux ~ 7.3 x 10~ n ergcm _2 s _1 ; and 
source A has a higher significance (24<r) with an energy flux of ~ 7.5 x 10 _n ergcm^ 2 s^ 1 
in the 0.2-300 GeV energy range (see Section PT7T1 the flux will be somewhat higher in 0.1- 
300 GeV). The use of different spectral models and different energy ranges in the two timing 
analyses may contribute to the discrepancy in the variability between source A and 3FGL 
J1638.6—4654. The spectrum of source A is fit to a power-law model and the spectrum of 
3FGL J1638.6—4654 is fit to a log-parabola model. Moreover, our timing analysis of source 
A uses photons in the energy range 0.2-300 GeV while 3FGL J1638.6—4654 is analyzed in 
the energy range 0.1-300 GeV. Our timing analysis would not be sensitive to flux variations 
(if any) below 0.2 GeV. 


3.3. Spatial Distribution Analysis 

We analyze the spatial distribution of source A, which is very likely to be associated with 
Kes 41, to examine whether it is a point-like or extended source. We apply both point-source 
and uniform-disk models with power-law spectra at the best-fit position to fit the emission 
in the energe range 2—300 GeV. In the point-source case, we set the spectral normalizations 
of the sources within 10° of Kes 41 as free parameters, and fix all the other parameters at the 
2FGL values. A TS value of 207 is obtained. In the disk case, the observed radius range for 
the uniform disks is 0?l-0?5 with a step of 0?1. We fix all of the spectral parameters of the 
sources at the values obtained above, but allow the spectral normalization parameters of the 


5 http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/aperture_photometry.html 
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disk models to be free parameters. The TS ex t value (calculated from 2 log(£disk/£pomt)) for 
each radius is smaller than zero, while the extended source detection threshold is TS e;r t = 16 
fjLande et al. 2012 ), which implies that no significant extended emission is detected. As a 
result, the GeV 7 -ray emission from source A seems to be point-like. 


3.4. Spectral Analysis 


The 7 -ray spectrum of source A is extracted via the maximum likelihood analysis of 
the LAT data in 6 divided energy bands from 0.2-300 GeV (see Tabled]). The spectral 
normalization parameters of the sources within 5° of Kes 41 are allowed to vary, but all of the 
other source parameters are fixed. In addition to the statistical uncertainties associated with 
the likelihood fits to the data, the uncertainty of the Galactic diffuse background intensity is 
considered. We vary the normalization of the Galactic backgroimd by ± 6 % from the best-fit 
values at each energy bin and estimate the flux from the object of inter est us ing the se new 
artificially frozen values of the background, following the treatment in Abdo et ahj (120091 ). 
The possible systematic errors are estimated to be 46% (0.2-0.5 GeV), 40% (0.5-1.0 GeV), 
20% (1.0-3.0 GeV), and < 15% (> 3 GeV). We keep only spectral flux points with TS higher 
than 4 (which corresponds to the detection significance of 2a) and derive 95% flux upper 
limits in the other energy bins. The obtained spectral data for source A are provided in 
Table dl 


We fit the 0.2-300 GeV spectral data of source A with a power-law model. The obtained 
spectral shape is relatively flat with a photon index of T = 2.38 ± 0.03. The flux is (9.2 ± 
1 . 0 ) x 10 -11 erg cm -2 s -1 , corresponding to a luminosity of ~ 1.6 x 10 36 d ^ 2 erg s -1 , where 
di 2 = d/12kpc is the distance to the MC asso ciated w ith SNR Kes 4 1 in units of the referecnce 
value estimated from the maser observation ( Koraleskv et al. 19981 . Also see Section |4~T1 for 
an estimate of the flux and luminosity with an exponential cutoff. 


4. DISCUSSION ON THE NATURE OF SOURCE A 

Based on our analysis of 5.6 years of Fermi-LAT data for the environment surrounding 
Kes 41, we have found a 7 -ray source detected at a significance of ~ 24a that appears to be 
coincident with the northwest rim of Kes 41. 

The relation between source A and Kes 41 is crucial for determining the origin of the 
7 -ray emission. In this section, we will discuss the possiblity of the 7 -ray emission arising 
from a pulsar and an SNR-MC hadronic interaction, respectively. 













4.1. A Pulsar? 


Galactic pulsars are important 7 -ray source can di dates, and there have been numerous 
pulsars detected by Fermz-LAT in recent years flAbdo et ah 2010e). Although the 3cr error 
circle here does not include any known pulsars, the possibility of correspondence to a pulsar 
associated with Kes 41 still cannot be ignored. Theoretically, there may be a descendent 
stellar compact remnant a fter t he core -collapse supernova (SN) explosion of the > 18M 0 
progenitor of the remnant flZhang et al.1120151 ). Such a compact stellar remnant has not been 
conclusively associated with Kes 41 in the literature. 

We fit the spectrum of source A with a power-law model with an expon ential cutoff, 
dAph/dAph = A'Ap h r exp(—Aph/A p h,cut), typical for a pulsar flAbdo et al . 2010a). The model 
fit yields A p h, C ut = 4.0 ± 0.9 GeV and the spectral index of T = 1.9 ± 0.1. In this model, the 
flux in energy range 0.2-300 GeV is (7.5 ± 0.9) x 10 -11 erg cm -2 s -1 , and the corresponding 
luminosity is (1.3±0.2) x 10 36 df 2 erg s -1 . The significance of the exponential cutoff power law 
(approximately described by yjTS cuto ffCr = JTSpj^r.utnff — TSpr ( cr) is ~ 6 ff. The^spectral 


shape of source A is similar to those of the detected 7 -ray pulsars flAbdo et al. 201 31). which 
usually show flat spectra below 1 GeV and exponential cutoffs in the energy range ~ 0.4- 
6 GeV. If this source is a “kicked” pulsar moving from the SNR center, then the best-fit 
position, 0.°1 away, would imply a projected traverse velocity of 180-4900di 2 kms _1 if the 
remnant’s age estimate 4-110 kyr flZhang et al.l 120151) is adopted. (The closer the position 
within the 3 a circle is to the SNR center, the lower the velocity would be.) The upper 
limit of the velocity seems very high, but there is also acuumulatiug evidence for high pul sar 
velocities, even exceeding 4 x 10 3 kms -1 (e.g., PSR B2011+38 and PSR B1718-35, Zou et al. 
20051 ). On the other hand, if it is an associated pulsar at ~ 12kpc, then its 7 -ray luminosity 


of the order of 10 36 er gs -1 (Sect ion 13. 4 j) would be among the highest among the (radio loud) 
7 -ray pulsars flAbdo et al.l 120131 ). which seems difficult to accept in view of the no detection 
of any radio pulsar here. 


4.2. Emission from Particles Accelerated by Kes 41? 


The 3cr error circle is on the northwestern boundary of SNR Kes 41 and essentially con¬ 
sistent with the shock-MC interaction region. Actually, it covers not only the 1720MHz OH 
maser but also the northwestern molecul ar gas at a systemic velocity of V LSR ~ —SOkrns ” 1 
that surrounds the remnant flZhang et al.l 20151 : see Fig. [3]). It is very possible that the 7 -ray 
emission arises from the relativistic particles accelerated by the SNR shock waves. We need 
to confront the leptonic and hardronic mechanisms with the obtained 7 -ray data. 
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4-2.1. Leptonic Scenario 


First, we consider the scenario in which the 7 -ray emission comes from the inverse 
Compton scattering off the relativistic electrons accelerated by the SNR shock. The emis- 
sivity of the bremsstrahlung process is compatible with that of the p-p process if the numbe r 
ratio of electrons to protons at a given energy, K ep , is of order ~ 0. 1 flGaisser et al.l 119981 ). 
Nevertheless, the values of K ep obs erved at Earth (I Yuan et ahl 2012) and predicted by the 
diffusive shock acceleration theory ([BelliIl978i ) are both of the order of ~ 0.01. Therefore, 
the bremsstrahlung gamma-ray is usually insignificant. 


We fit a power-law electron spectrum with a cutoff, dN e /dE e oc E~ ae exp (—E e /E e:CUt ), to 
the spectral data and only consider the cosmic microwave background as the seed photons 
(referred to as Case A). As can be seen in Fig. [4] (blue dotted line), the fitting effect is 
less satisfactory. We obtain a e ~ 2.0 and E Sjcu t ~ 400 GeV (also see Table [2]). The 
normalization is given by the total energy deposited in electrons with energy above 1 GeV, 
W e (> lGeV) ~ 1.3 x 10 51 erg. This electron energy budget is unreasonably high as the order 
of the canonical SN explosion energy. 


4-2.2. Hadronic Scenario 

Next, we consider the scenario in which the 7 -ray emission is produced by the collision 
of the shock accelerated protons with dense molecular gas. For the case (referred to as Case 
B) in which the protons collide with the dense target molecular gas (with average number 
density n t ), we assume for the protons a broken power-law distribution, dN p /dE p oc E p “ p (l + 
(E p / £’b) 2 )~ A “ p//2 , to fit the spectral data (see Fig.[4](red dashed line) and Table[2|). We thus 
obtain a p = 2.0, Aa p = 1.2, and a break energy of E h = 18 GeV. The total energy deposited 
in the protons with energy above 1 GeV is W p (> 1 GeV) ~ 0.7 x 10 50 iAi(nt/100 cm -3 ) -1 
erg, where E 5 1 = F^n/IO 51 erg is the dimensionless SN explosion energy. SNR Kes 41 
has been found to be surrounded by molecular gas of density with uh 2 ~ 140-500 cm -3 
in the northwest and HI gas of density with n(HI) ~ 40 cm -3 in the southeast (also see 
Fig. |3]). If the mean target density n t is approximately of the order of 100 cm 3 , then 
W p ( > 1 GeV) ~ 1 x 10 50 £ 5 i erg, namely, the fraction, 7 , of the SN explosion energy converted 
to protons is of the typical order of 0.1. While in this scenario the hadronic 7 -rays are emitted 
at the SNR shock, it is noteworthy that the centroid of the 3<r circle of source A appears to 
be outside the northwestern boundary of the SNR. 

The hadronic emission can alternatively be considered as originating from the adjacent 
MCs that are “illuminated” by the diffusive relativistic protons escaping from the SNR shock 
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front. In the finite volume of a nearby cloud, the protons’ energy distribution can be obtained 
by calculating the diffusive escaping protons accumulatively throughout the his tory of the 
SNR expansion (Case C ). For such a calculation, in the following, we refer to Li & Chen 
(.2012) and the references therein for details of the model. 


In the model calculation, we assume a converted CR proton energy fraction of rj = 0.1 
and an SN explosion energy of E SN = 10 51 erg. The SNR radius in the southeast-northwest 
orientation is adopted as R s ~ 11 pc. The y-rays are assumed to arise from an MC, of 
thickness A R c , which is in contact with the shock surface; the refore, t he MC c enter is at 
R c = R s + AR c /2 away from the SNR center. According to Zhang et al. ( 2015 1. the SNR 
evolves in a cavity and may have been drastically decelerated and entered the radiative 
p hase as soo n as the blast wave encountered the cavity wall, after a Sedov evolution lifetime 
( Sedov 1959 1 of t enc = 4 x 10 3 (n H /0.3cm _3 ) 1/,2 £ , 5 ’ 1 1 ^ 2 (R s /ll pc) 5 / 2 years. We assume that 
the particle acceleration process is not significant after this time. Therefore, the average 
distribution of the cumulative escaping protons in the volume of the MC at the remnant age 
t age is rewritten as 

F ave (E p ,t age ) = / r 2 dr / / / f(E p , R het {R c ,ti, 9, (ft), tax) (t;)sin6> dd dcftdti 

jRc-ARc/2 Jo Jo Jo 

rR c +AR c /2 

/ / r 2 dr , (1) 

' Jr c -ar c /2 

where t* is the time at which a proton escapes from the SNR shock, Cif = Cge — C is the 
diffusion time after escape, Abet is the distance between the escape point on the shock surface 
and a given point in the cloud (with position angles (9, (ft)), and f (E p , R het (R c ,ti, 6, (ft),t dd ) 
is the distribution function at a given point of the protons that escape from the unit area 
at an arbitrary escape point. Considering the remnan t’s age range ^4- 100 kyr estimated 
from the ionization timescale of the X-ray emitting gas (lZhang et al.ll2015n . we calculate the 
model with three age numbers, 4, 10, and 100 kyr. This model can fit the spectral points as 
well, as exemplified by the solid line for i age = 10 kyr in Figure SI The model parameters are 
listed in Table [21 The photon index a p = 2.4, the energy-dependent index of the diffusion 
coefficient 6 = 0.7, and the correction factor of slow diffusion around the SNR \ ~ 0.01-0.1 
are in normal ranges. The A R c value ~ 5-13 pc (~ 0?02-0?06) implies that the MC involed 
in the p-p hadronic interaction is essentially within the 3<r circle of source A. Note that such 
a source size is much smaller than the PSF size 0?5 of the Fermi- LAT at energies above 
2GeV, consistent with the above judgement of a point-like source. The “illuminated” MC 
mass M c i, around 1O 5 M 0 , seems reasonable as compared with the mass of the molecular gas 
“reservoir” in the northwest, which is no less than ~ a few times 10 5 M ^, which is estimated 
from a limited field of view of the CO observation (lZhang et al.l 1201 51). 
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However, th e cavi ty wall may send a reflected shock backward when the blast wave col¬ 
lides with it flZhang et ah 2015). If the reflected shock can still effectively accelerate particles 
after the forward shock becomes radiative, then the situation would be more complicated 
than the above cases. For simplicit y, we approx imate this case as a continuous proton in¬ 
jection from the SNR center (lAharonian fe Atovan 19961 ) (Case D ). In this case, the MC 
is regarded as a point at R c from the SNR center and the same energy conver sion fraction 
r] = 0.1 is adopted. We follow the algorithm described in Aharonian fe Atovani (1l996n and 
fit the spectral data, as exemplified by the green dashed line for t age = lOkyr and R c = 20 pc 
in Figure [4] These model results are generally similar to those of Case C, with a slightly 
harder model spectrum at > 100 GeV. For the three sets of parameters with f age = 10 and 
100 kyr, we again have a p = 2.4 and S = 0.7. The x values are ~ 0.05-0.5 in a normal 
range. A higher mass of the “illuminated” part of MC than Case C is required, but is still 
consistent with the MC mass estimate from the CO observation in the order of magnitude. 


The hadronic scenarios, both the interaction at the shock (Case B) and the illumination 
by escaping protons (Case C/D), can generally explain the y-ray properties of source A. The 
escape cases have harder model spectra at > 10 GeV than the interaction-at-the-shock case. 
Further TeV observations will likely be of help to distinguish the two scenarios. 


5. COMPARISON WITH OTHER GeV-DETECTED SNRs IN MC 

ENVIRONMENTS 


We now present a brief discussion of Kes 41 within the context of other Galactic SNRs 
that have been detected at y-ray energies. An intriguing trend has emerged in these studies 
where the Galactic SNRs that are known to be interacting with dense clouds and that are 
detected at (very) high energies also appear to exhibit contrasting morphologies in the X-ray 
and the radio. Specifically, these sources exhibit the shell-like radio morphologies that are 
characteristic of SNRs coupled with a center-filled X-ray morphology that is thermal in origin, 
and therefore belong to the class of thermal composite or mixed-morphology SNRs (also see 
Section |T]) . Actually, about half of the 36-37 known thermal com posit es have been found to 
be interacting with adjacent MCs (see Table 4 in Zhang ef ah 20151) . While the origin of 
these contrasting morphologies remains uncertain, it appears that the interaction between 
the SNRs and the dense clouds plays a crucial role. Proposed origins for these morphologies 
include the evaporation of shock-engulfed cloudlets, thermal conduction within the interior 
h ot gas, and heati ng by the shock reflected from the wind-cavity wall; the reader is referred to 
Chen et al.l (120081 ) and references therein for a detailed review of these proposed mechanisms. 


We tabulate the thermal composite SNRs that have been detected at y-ray energies by 
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Fermi -LAT in Table [3j So far, there are 13 SNRs (including Kes 41) of this class that have 
associated GeV 7 -ray emission, and an additional six of them possibly have associated GeV 
7 -ray emission, as listed in Table [3J We can see that most of the GeV-detected thermal 
composites are in physical interaction with MCs. 


In Table |3l we collect the photon indices in the GeV band and adopt the luminosities in, 
or convert them to, the 1-100 GeV energy range for ease of comparison. The ~GeV spectra 
of these SNRs are soft, with power-law photon indices of T > 2.0, in distinct contrast with 
the hard spectra (T ~ 1.4-1. 8 ) of the supp osed leptonic process dom i nated 7 -ray SNRs, 
e.g., RX J0852.0-4622 ( Tanaka et al. 2011 1 and RCW 86 ( Yuan et al. 2014). Except for 
HB 21 and Kes 27, the 1-100 GeV luminosities of the 13 identified GeV 7 -ray sources are 
on the order of a few times 10 35 ergs -1 , which are significantly higher than t hose o f the 

< 10 34 ergs -1 for RX J0852.0-4622, Tanaka et ah 


2011 

; and RCW 86 , 

Yuan et al. 

2014) 


Yuan e_t ah 201411 . For exceptional cases of HB 21 and Kes 27, the low 


lumi nosities may be due to proton collisio ns with only a very small amount of dense clouds 
(e.g., IPivato et al.l 12013c King et al.l 1201511 . We note that where detailed modeling has been 
applied to the 7 -ray spectra of these sources, hadronic models have generally proved to give 
better fits to the data than leptonic models (except for the uncertain cases of Kes 17 and 
HB 9). 


These past 7 -ray observations of thermal composite SNRs-including the observation of 
Kes 41 that is presented in this paper - have thus produced insights into how SNRs interact 
with MC and how SNRs accelerate CR particles. Additional 7 -ray observations of thermal 
composites are necessary and timely to explore the relation between emission at these high 
energies and the origin of the contrasting morphologies that characterize SNRs of this type. 


6. Summary 

We perform an analysis of the 7 -ray emission in a 14° x 14° region centered on the 
thermal composite SNR Kes 41, using 5.6 years of Fermz-LAT observation data. We End 
a point-like source to the northwest of the SNR with a significance of 24<j in 0.2-300 GeV. 
Neither significant long-term variability nor periodicity is detected from the timing analysis 
of source A in the same energy range. The 3er error circle, 0?09 in radius, covers the 1720MHz 
OH maser and is essentially consistent with the location of the V LSR -50 kms -1 MC with 
which the SNR interacts. The source emission can be described by a power-law spectrum with 
an exponential cutoff with a photon index of 1.9 ± 0.1 and a cutoff energy of 4.0 ± 0.9 GeV. 
The corresponding 0.2-300 GeV flux is (7.5 ± 0.9) x 10 - 11 ergcm - 2 s -1 , and the luminosity 
is ~ 1.3 x 10 36 ergs -1 at a distance of 12 kpc. Although the spectrum is similar to those of 
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pulsars, there is no radio pulsar in the 3a circle responsible for the high luminosity. While 
the power-law electron spectrum with a cutoff for inverse Compton scattering would lead to 
a difficulty in the electron energy budget, the emission can be naturally explained by the 
hadronic interaction between the relativistic protons accelerated by the shock of SNR Kes 41 
and the adjacent northwestern MC. By comparison with the hadronic interaction at the 
shock, which appears off the best-fit position of the source, illumination of the adjacent MC 
by the protons escaping from the shock front seems more consistent with observations. A 
list of Galactic thermal composite SNRs detected at GeV y-ray energies by Fermi -LAT is 
presented in this paper. 

B.L. is grateful to Xia Fang, Ning-Xiao Zhang, and Zheng-Gao Xiong for the help about 
Fermi data analysis. We thank the support of NSFC grants 11233001 and 11403075. This 
work has also benefited from 973 Program grant 2015CB857100, grant 20120091110048 from 
the Educational Ministry of China, and the grants from the 985 Project of NJU and the 
Advanced Discipline Construction Project of Jiangsu Province. This research has made use 
of the SIMBAD database, operated at CDS, Strasbourg, France. 
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Table 1. Fermi LAT Flux Measurements of Source A in the Kes 41 Region 


E ph (energy band) 
(GeV) 

E* h dN(E ph )/dE ph a 
(10- 12 erg cm -2 s 1 ) 

TS value 

0.32 (0.20-0.50) 

15.9±3.3±7.3 

39 

0.71 (0.50-1.00) 

24.4±4.6±9.9 

113 

1.73 (1.00-3.00) 

21.3±1.8±4.6 

224 

5.48 (3.00-10.0) 

10.7±1.2±1.7 

105 

17.3 (10.0-30.0) 

2.4±1.0±0.3 

9 

94.9 (30.0-300) 

< 2.9 b 

2 


a The first column of errors lists statistical errors and 
the second lists systematic errors. 

b The 95% upper limit. 


Table 2. Model Parameters for the Emissions of SNR Accelerated Particles 





















- 18 - 


Table 3. 


Parameters of the 7 -ray Emission of the Galactic Thermal Composite SNRs 


Obtained from Fermi -.LAT Observation 


Source 

Distance 

r 

^1-100 GeV 

MC interaction 

References 


(kpc) 


(10 35 ergs -1 ) 




G6.4-0.1(W28) 

2.0 

2.74 ± 0.06 b 

1.0 

Y 

1,2 

G31.9+0.0(3C 391) 

7.2 

2.50 ± 0.04 b 

4.0 

Y 

3,4 

G34.7-0.4(W44) 

2.8 

3.02 ± 0.10 b 

2.7 

Y 

5,6 

G43.3-0.2(W49B) 

8 

2.29 ± 0.02 c 

8.0 

Y 

7,8 

G49.2-0.7(W51C) 

6 

2.5 ± 0.1 b 

4.4 

Y 

9,10 

G89.0+4.7(HB 21) 

1.7 

2.33 ± 0.03 c 

0.13 

Y 

11,12 

G189.1+3.0(IC 443) 

1.5 

2.61 ± 0.04 b 

1.0 

Y 

5,13 

G304.6+0.1(Kes 17) 

9.7 

2.0 ± 0.3 C 

12 

Y 

14,15 

G327.4+0.4(Kes 27) 

4.3 

2.5 ± 0.1 c 

0.24 


16,17 

G337.8-0.1(Kes 41) 

12 

2.38 ± 0.03 c 

7.7 

Y 

18 

G348.5+0.1(CTB 37A) 

11.3 

2.19 ± 0.07 c 

7.8 

Y 

19,20,21 

G357.7-0.1(MSH 17-39) 

12 

2.5 ± 0.3 C 

5.8 

Y 

22,23 

G359.1-0.5 

7.6 

2.60 ±0.05° 

4.0 

Y 

24,25 


G0.0+0.0(Sgr A East) (?) d 

8.0 

2.32 ± 0.03 c 

8.7 

Y 

26,23 

G132.7+1.3(HB 3) (?) d 

2.2 

2.30 ± 0.11 c 

0.04 

Y? 

27,23 

G156.2+5.7 (?) d 

3 

2.35 ± 0.09 c 

0.26 


28,23 

G290.1-0.8(MSH 11-61A) (?) d 

7 

~ 2.28 c 

1.5 

? 

22,23 

G160.9+2.6(HB 9) (?) e 

1.0 

2.30 ± 0.05 c 

0.013 

? 

29 

G166.0+4.3 (?) e 

4.5 

2.27 ± 0.1 c 

0.11 

? 

30,31 


Note. — (1) Frail [201 lJ: (2) Abdo et a 

. 2010al: (3) iRadhakrishnan et al. 19721: (4) Ergin et al. 2014 

; (51 

Seta et al.lIl998l: (6llAbdo et al.ll2010d; (7) 

Brogan Sz Troland 200ll; (8) lAbdo et al. 2010bl: 

(9) iKoo et al.l 

2005 

(10) 

Abdo et al. 2009|: (ID Byun et alj 20061: (12) Pivato et al. 

2013 (131 Abdo et al. 

2010d 

; (14) 
(18 

Combi et al. 

201C 

: (lsicielfand et al. 20131: (16lMcClure-Griffiths et al.l 

200lL (17lXing et al. 

12015 

Zhang et al. 

2015 

(19lRevnoso & Manguml l2000l: (20' 

Oastro & Slanel 201C 

; (2i' 

Brandt & Fermi-LAT Collaboration 20131: 

(22 Rosado et al.1 19961: (23lAcero et al. 

20151: (24lUchida et al. 

19921; (25lHui et al 

1 I 2011 

; (26lReid 

19931; (27 iRoutledge et al. 199ll; (28lReich et alJ 19921: (29lAraval 20141:(30'Landecker et al. 19891: (31 

Aravall20li: 


a Adopted from [jiang et al* l l2010l SNR-MC association table. 


b The photon index above the break energy for broken power-law spectrum. 
c The photon index of single power-law spectrum. 

d Question mark: the association of the detected 7-ray emission with the SNR is uncertain. 
e Question mark: not listed in the latest 3FGL catalogue llAcero et al.|[2015h . 
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Fig. 1.— TS map (2-300 GeV) of 1° x 1° region centered at Kes 41. All sources except 
2FGL J1638.0—4703c have been subtracted. The green cross labels the position of a 2FGL 
source, the cyan cross labels the position of 2FGL J1638.0—4703c, and the dashed blue circle 
indicates the 3cr error range of the best-fit position for the residual emission found in the 
Kes 41 region. The image is overlaid with the MOST 843 MHz radio c ontou rs (in green) (at 


seven linear scale levels between 0.00 and 0.79 Jy beam 1 ; from IWhiteoak fe Green 


The magenta diamond represents the OH (1720 MHz) maser spot (IKoraleskv et al. 


(19961 


19981 ). 
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Fig. 2.— Monthly 7 -ray light curve of source A in the energy range of 0.2—300 GeV 
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Fig. 3.— Tri-color image of Kes 41 in multiwavelengths. Red: Fermz-LAT 2-300GeV 
counts map centered at SNR Kes 41, smoothed with a Gaussian of width 0?6 (per pixel bin 
representing 0?01). Blue: 12 C0 (<7=1-0) integrated emission (Vlsr = — 70 to —40kms _1 ) 
with a field of view of 11' x 10'. Green: H I line emission from SGPS integrated map 
(Ksr — —55 to —50 km s' 1 ). The green contours and the green and cyan crosses are the 
same as in FigjJ] The white curves show the TS = 100,144 and 196 contours (which 
correspond to significance 10cr, 12<r and 14cr, respective ly). The white diamond indicates 
the location of the 1720 MHz OH maser flKoraleskv et al.1119981 ) and the green circles label 
the positions of known pulsars. The dashed blue circle indicates the 3cr error circle of the 
best-fit position of source A. 
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Fig. 4.— Fermi 7 -ray spectral energy distribution of source A fit with various models (see 
text). Systematic errors (see Section 13.4(1 are indicated by black bars and the statistical 
errors are indicated by red bars. 

















